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' Abstract. Recent developments in laser physics have called renewed attention to Laguerre- 

Gaussian (LG) beams of paraxial light. In this paper we consider the corresponding 
LG modes for the two-dimensional harmonic oscillator, which appear in the transversal 
plane at the laser beam's waist. We see how they arise as Wigner transforms of Hcrmite- 
Gaussian modes, and we proceed to find a closed form for their own Wigner transforms, 
m : providing an alternative to the methods of Simon and Agarwal. Our main observation 

is that the Wigner transform intertwines the creation and annihilation operators for the 
two classes of modes. 
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1. Introduction 

■ 

■ In the early 1990s it was observed that a Laguerre-Gaussian (LG) beam of paraxial 

light has a well-defined orbital angular momentum and that such a beam may be created 
from a Hermite- Gaussian (HG) beam by way of an astigmatic optical system pQ. This 
discovery was soon to find applications, for example, in biology (for "optical tweezers") 
and in the study of quantum entanglement [2J. Additionally, recent experiments suggest 
that beams with orbital angular momentum might be generated in situ, in free electron 
q | lasers j3]. 

In this paper we consider the corresponding LG modes for the two-dimensional harmonic 
oscillator, which appear in the transversal plane at the laser beam's waist. Indeed, one 
may use the operator algebra of the harmonic oscillator to study the analytical forms of 
the HG and LG beams [1], [5]. Here we demonstrate that the LG modes arise as Wigner 
transforms of HG modes, and we proceed to find their own Wigner transforms. Our 
methods provide an alternative to those of Gase [B] and those of Simon and Agarwal [7]; 
in particular, Simon and Agarwal were first to discover the closed form expression that 
we rederive here, using our new point of view. 
We define the d- dimensional Wigner transform as 



X 



for functions /, g of d variables. If d = 1, we omit the subscript. The Wigner transform 
is often restricted to the case when / = g, and then one writes Wd(f) '■= Wd(f, /). More 
generally, we define the extended Wigner transform, of a function F of 2d variables, as 



W d (FM£) = (2k)-* J e**F 
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This more general transform is fundamental for the study of LG modes; this is our main 
insight and provides the central theme for our paper. 

Our main results are summarized in the following theorem. 

Theorem 1. Let hj denote the j th Hermite function, so that hjk(x,y) = hj(x)hk(y) are 
the HG modes. 

(a) The extended Wigner transform intertwines the creation and annihila- 
tion operators of the HG and LG modes. 

(b) The LG modes are precisely W(hjk): the extended Wigner transforms of 
the HG modes. 

(c) The Wigner transforms of the LG modes are given by 
W 2 (W(h jk ),W{h mn ))(x,i) 



W {hjm ) — j W(K„) 



We remark that part (b) of the theorem follows from part (a) and the fact that W(hoo) = 
hoo- 

In the case (j, k) = (m,n), the formula ([I]) was proven by Simon and Agarwal, us- 
ing the metaplectic representation, and in fact their proof extends to the general case 
without modification [7J. The generalization should be useful in light of the relationship 
between the Wigner transform and the Weyl quantization of observables. Indeed, for 
any tempered distribution a e iS'(IR 2d ) we may define the Weyl quantization of a as the 
[pseudo differential] operator given by 

(Op»«) (x) = 2---- d // exp {^^-) - (^*) <V) 

The normalization factor is chosen so that Op H/ (l) is the identity operator. Then it is 
easy to check that for any Schwartz functions /, g G 5(IR d ) we have 

(f\0 P w (a)g) = 2- d n-i J J a(x, £)W d (f,g)(x, dxd^. 

For an exposition of the Weyl calculus of pseudodifferential operators, the reader may 
consult the book of Gerald Folland [8]. 

In Section [2] we review the theory of the one-dimensional harmonic oscillator. Since 
this subject is very well known, we will be brief, so the section is mostly a means of 
setting up notation for the following sections. In Section [31 we review the theory of the 
two-dimensional harmonic oscillator. Our primary reference for Sections [2] and [3] is the 
classic book of Messiah [9]. 

In Section H] we write the {(^-representation in terms of complex variables, which 
provides a clean expression for the LG modes in terms of Laguerre polynomials. The 
connection with LG modes of paraxial light is briefly exhibited in Section 

The LG modes are expressed as Wigner transforms in Section [HI and in Section [7] we 
calculate their own Wigner transforms. 
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2. The One-dimensional Harmonic Oscillator 
We begin by setting up the notation for the one-dimensional harmonic oscillator, given 

1 



We set 



so that 



H = ^(P 2 + Q 2 ), where [Q,P] 



a=-^(Q + iP) and a f = -^=(Q - iP) 



V2' 



1 



Moreover, we define 
so that 



[a, at] = 1 and H = -(aa) + a) a). 



N = a) a 



H = N + 



It is well known that Spec(iV) = {0,1,2,...} consists of nondegenerate eigenvalues 
and that at and a are the "ladder operators": the creation and annihilation operators, 
respectively. The associated eigenvectors form a complete set, so we can normalize to get 
an orthonormal basis of eigenvectors for the observable N: 

|0>, |1>, |2>,... 

corresponding to the eigenvalues 

0, 1, 2,... 

In the {Qj-Representation (that is, the Schrodinger representation where states are 
considered as functions of position) the ladder operators are written as 



a = 



V2 



x 



and 



The ground state ho satisfies 



so that 



1 / 



dx 
d 



y/2\ dx) ' 



d 

dx 



+ X 



h {x) = 0, 



h (x) = it *e 2X . 

Using the ladder operators, we get the rest of the eigenvectors, which are precisely the 
Hermite functions: 



h n (x) = it 2 2 2 l x 



d 

dx 



= 7r-3(n!)-2 2 -f(-l) n e^ 



e 2 



d n 
dx n> 



it *(n\) 22 2 6 2 X HJx) 
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where the H n are the Hermite polynomials 



,2 d n _ x 2 



H n (x) = {-We* 



3. The Two-Dimensional Isotropic Harmonic Oscillator 
In two dimensions, we have 

H= l -{Pl + Pl + Q\ + Ql) 
with creation and annihilation operators inherited from the one-dimensional case: 

aj = ^(Qj + iPj), J = 1,2. 

These are interpreted as creation and annihilation operators, respectively, of quanta of 
type j E {1,2}. They satisfy 

[a,, at] =Sij. 

And we have the corresponding "number operators" 

N j = a]a j , j = 1,2, 

with 



representing the total number. 
Now the eigenvalues for 



N = Ni + N 2 



H = N 



are 



{1,2,3,...} 

and the eigenvalue j has degeneracy j. Hence there are many possible choices of eigenbasis. 
One often simply takes the basis consisting of tensor products of the one- dimensional 
eigenvectors, corresponding to the complete set of commuting observables {N 1: N 2 }. The 
elements of this basis are given by 

\ ni n 2 ) = (n 1 !n 2 !)-^a t 1 ni 4 n2 |00). 

These are precisely the HG modes; in the Schrodinger representation, they are tensor 
products of Hermite functions. 

But here we will construct another basis. 

The angular momentum operator L is defined by 

L = QiP 2 - Q-2.Pi = i{a-ya\ - a[a 2 ), 

and one may check that it is a constant of motion. We will show that {A^, L} is another 
complete set of commuting observables. 
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Let 

A± = -j={a 1 Tia 2 ) 

and 
Then 

[A r ,A s ] = [Al,Al] = and 
[A r ,Al\=5 rs for r,se {+,-}. 

We think of A\, A r as creation and annihilation operators, respectively, of quanta of 
type r e {+, — }, so then 

N r = A\A r 

represents the number of "r quanta" . 

Hence the problem of finding eigenvectors common to iV + and iV_ is formally equivalent 
to finding eigenvectors common to Ni and N 2 . So by the usual arguments we see that 

Spec(7V+) = Spec(AL) = {0, 1,2,.. .} 

and that these two observables, N + and N_, form a complete set of commuting observ- 
ables: to each (n + , there is a common eigenvector, denoted by that is unique 

to within a constant. 
In fact, 

A + \00) = A_\00) = 

and the states 

|n + n_) = {n + \nJ)-^A ] l l+ A ] l l -\m) 
form a complete orthonormal eigenbasis common to iV + and iV_: 

N + \n + n_) = n + \n + nJ) 

and 

N^\n + n_) = n^\n + n_). 

We find that 

N = N + + N_ 

and that 

L = N + -N_. 

Hence and L form a complete set of commuting observables. Moreover, 

[L, A{] = ±4 

and 

[L,A±] =TA±, 

so when they act upon an eigenvector of L, A\_ and A_ increase L by one unit, and A_ 
and A + decrease L by one unit. So it is natural to consider iV + as the number of particles 
with positive charge, iV_ as the number of particles with negative charge, and L as the 
total charge (to within a constant). 
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4. The Two- Dimensional Isotropic Harmonic Oscillator in the 

{(^-Representation 

In the {(^-representation, we write 



(2) 



so that we have 



,t_ w, d \ j 1 f 9 



(3) 



1 / 9 . / <9 

i / 9 / d 

A -=2{ X+ d- X +l { y+ d-y 

a\ 1 ( 9 f 9 

A+= 2\ X -d- X + \ y -d-y 



Al 



If d . ( d 



2 \ X dx 1 y dy 



In complex notation, these take a simple form. We write z = x + iy, so that 

d _ 1 / d 1 d \ d _1 / d _1 d 

dz 2 \dx i dy J &z 2 \ dx i dy 



Then 



9 1 9 



+ ~2 5z' -"2 0*' 



Moreover, we can write 



; 9 , .4. 3 _1. 

'— e 2 and = -e2 zz — e 2< 

oz dz 



The ground state is given by 

u (z,z) = -K~^e~^ z * = (z,z\00), 

and of course we get all other eigenvectors by applying the creation operators to this. 
These are precisely the LG modes. Explicitly, 

(z,z\n + n_) = (z,z\{n + \nJ^A^ + A^-\m) 

= ,-W-o-»(-ir++"-^(|)" + (£) 

We will now formulate this in terms of Laguerre polynomials. 



n+ / o \ fi- 
zz 



e 
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First suppose that n + >n_. Then 

(z,z\n+n-) = 7r _ s(n + !n_!)-3(-l) n++n -e3 ilz f-^j [(-z) n +e- zT ] 

= 7r-5( n+ ! n _!)-|(-l)«- e 2^(™— »+) (jr^—) [(zz) n+ e- z "] . 



d{zz) 



We define the Laguerre polynomials by 



nV 7 n! dx n v 7 
for n > and a > — 1. Then we see that, when n+ > n_, 



i / ri-i 



+ 

i 
■i 



,z,z\n + n_) = 7T-5 (_i)«-2K-"-) e 4"L«+-»-(^). 
On the other hand, if n + < n_, similar arguments show that 

(^|n + n_> = tt"* f^±j) 2 (-l)"+^-"+)e-^L--"+(^). 
We summarize this in the following theorem: 

Theorem 2. We consider the states in the {Q} -representation with position co- 

ordinates (x,y) G M. 2 , and we write z = x + iy. Then 



(z, z\n + n_) 




(-l) n -z (n +- n ^e-^ zz L™+- n ~(zz) ifn + >n_, and 



5. Comparison With LG Modes of Paraxial Light 

The time-harmonic field amplitudes of LG optical modes are expressed in cylindrical 
coordinates, with propagation coordinate z along the beam axis, as 



itp/r) oc e- i${r) e ^ 




where the phase is 



$(r) = icj) - kz + — — -(2p + £+ l)tt(z) 
2R{z) 



and where again l)p is the Laguerre polynomial with radial and azimuthal indices p and 
£, respectively. The term W(z) = tan _1 (^) is part of the Gouy phase, R(z) = Zr ~^ z 

is the radius of curvature, w(z) = Wq^Ji + (j^) 2 is the beam waist, and z R = ^kw^ is 

the Rayleigh range. We see that the total phase evolves helically along the propagation 
coordinate z. For more on this, see, for example, [2], [3], [10J. 
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In the transversal plane at the beam waist (z=0), we get 

which is essentially what we have in Theorem [2j 

6. LG Modes Are Wigner Transforms 

Physicists in the optics community recently computed the Wigner transforms of LG 
modes [6], [7]. Surprisingly, the states |n+n_) are already themselves Wigner transforms 
of Hermite functions. The Hermite functions, as in Section [21 are given by 

hjix) = 7r-i(j!)-§2-i(-l)a* 2 |U-* 2 , 

which, as we know, satisfy the orthogonality relations 

(hj\h k ) = Sjk- 
And we continue to write the HG modes as 

hjk(x,y) = hj(x)hk(y). 
Now, for reference, we list the standard properties of the Wigner transform: 



(4) W(f,g) = W(g,f), 

(5) j W(f,g)(x,Od£=V^f(^g^y 

(6) J J W(f,g)(x,0dxd£ = 2V^(f\g), 

(7) jw(f,g)(x,Odx= y/tof (J^=\g (J=) , 

where /(£) = ^7= j e~ lx ^f(x) dx is the Fourier transform, and ("Moyal's identity") 

(8) (W{h, gi )\W(f 2 ,g 2 )} = JMN)(gi\g 2 ). 



These properties are well known and are easily checked. But there is another elementary 
property which seems to be a new observation: 

Theorem 3. The extended Wigner transform intertwines the creation operators of the HG 
and LG modes, and it intertwines the annihilation operators of the HG and LG modes: 

A ] + W = Wa\ , AlW = Wa\ , 

A + W = Wa u A.W = Wa 2 . 
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The proof is a direct calculation, using the expressions (jSJ) and (j3J). As a first application, 
it is now easy to see that the extended Wigner transform commutes with the operators 
H and N. 

Additionally using the fact that W(h 0Q ) = h 00 , we have the following corollary: 

Corollary 1. We consider the LG modes |n + n_) in the {Q} -representation with position 
coordinates (x,y) £ M 2 . Then we have 

(x,y\n + n^) = W(h n+ , h n _)(x,y). 
Combining this with Theorem [2J we have a formula: 
Corollary 2. Suppose x, y £ R, and let z = x + iy. Then 

W(hj,h k )(x,y) = 



i 

7T~5 (Mj 5 (-l) fc ^'- fc e -5^L^ fc (^) z/ j > k, and 



Corollary [2] by itself is not new. For example, two independent proofs may be found 
in Folland's book [S], although we caution the reader that his statement of the result 
contains a misprint. 

So we have yet another way of understanding the orthogonality relations of the states 
\n + nJ): as being inherited from the orthogonality of the Hermite functions. Explicitly, 

(m + m_|n + n_) = (W(h m+ , h m _)\W(h n+1 h n _)) 
(hm+ 1 h n+ ) I h n _ ) 

As another application of Corollary [Tj, we can think of the LG modes \n + n_) as being 
"interference terms" resulting from the quadratic nature of the Wigner transform. If 
we apply the Wigner transform to sums of Hermite functions, we get interference terms 
that are LG modes with nonzero angular momentum. In fact, we have the following 
polarization identity: 

(x,y\n + n-) = ~W(h n+ + h n _)(x,y) - ~W(h n+ - h n _)(x,y) 

+ -^W(h n+ - ih n _)(x, y) - - W(h n+ + ih n _)(x, y). 

One may think of the Wigner transform as being a kind of "product". Just as Hermite 
functions give rise to two-dimensional HG modes, via the tensor product, Hermite func- 
tions also give rise to two-dimensional LG modes, via the Wigner transform. To complete 
this series of relationships, Simon and Agarwal [7] noted that the two-dimensional HG 
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and LG modes are unitarily related by the operator 

exp —T\ , where 2\ = xy — 



4 / dxdy 

Alternatively, we can use the Wigner transform to again show that the HG and LG 
modes are unitarily related. For functions F of two variables, we consider the extended 
Wigner transform of F: 

If we write the ^-rotation operator as 

(KlF)(x,p) = F(R«(x,p)) 

' x — p x + p" 

and if we write the partial Fourier transform as 

{F 2 F){x,y) = -L / e~^F(x,p)d Pl 

then we may write the extended Wigner transform as 

(9) WF(x,y)=(F 2 m 4 F)(x,y). 

This is clearly a unitary operator, so, in particular, we again see that the HG modes hjj. 
are unitarily related to the LG modes W(hjk). 

The results of this section are surprising, despite their mathematical simplicity. We are 
accustomed to thinking of the Wigner transform of a function as being a phase space rep- 
resentation of the function, but here the relevant Wigner transform is actually a function 
of the position vector (x, y) G M 2 , as it is an LG mode. So if we use the Wigner transform 
to study the phase space properties of LG modes, we are really taking Wigner transforms 
of Wigner transforms. Gase [BJ, followed by Simon and Agarwal [7], already recently 
derived expressions for Wigner transforms of LG modes (Gase in terms of a quadruple 
sum, and Simon and Agarwal in terms of a closed form involving no summation at all). 
We will do the same in the next section, but now using our new point of view. 



7. Wigner Transforms of LG Modes 

We now proceed to compute the Wigner transforms of the LG modes W(hjk)- In light 
of ([H]), the first step is to consider the action of the Wigner transform on partial Fourier 
transforms. 



Lemma 1. 
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This lemma follows from a direct calculation. 

We now apply the lemma to 

F = Rtt_ hjk and G = HX h mn . 

4 ■> 4 

That is, 

vf-\ u f x i~ X A ? (x 1 + x 2 

F(l) = H~^rv ** wr 

and similarly for G. Hence the Wigner transforms of the LG modes are given by 
W 2 (W(h jk ),W(h mn ))(x,0 

1 f f e *s-<** h . ( (xi + s) + (6 + *) ^ ^ ( (xi + s) - (6 + t) 



2« J J J \ 2 J \ 2 



We may simplify this integral by again rotating the variables by |. We state the result 
as a theorem. 

Theorem 4. The Wigner transforms of the LG modes are products of LG modes, given 
by the formula 

W 2 (W(h jk ),W(h mn ))(x,0 

= W{h jm ) -=—, -=— W(hkn 



In the special case (j, k) = (m, n), we have the following formula of Simon and Agarwal 
[7j . It is now proven by combining Corollary [2] with Theorem HI 



Corollary 3. For = {x x , x 2 , £ 2 ) G M 2 x M 2 ; we Zei 

1 
2 



TTien we /iai>e 

w^ 2 (^(/ ljfc ), w(/i ifc ))(x, 15 = Tr^-iy+v^L^Qo + g 2 )^(g - Q2) 

For comparison, we also give the Wigner transforms of HG modes: 
W 2 (h jk ,h mn )(x,i) = W(h jm )(x 1 ,^ 1 )W(h kn )(x 2 ,^ 2 ). 
In the particular case (j, k) = (m, n), if we let 

Q, = \{x\-xi + e 2 -? 2 2 ), 

then we have the formula 

W 2 (h jk , h jk )(x,0 = n-\-iy +k e- Q °L%Q + Q 3 )L° k (Q - Q 3 ). 
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